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Abstract 

We apply Pauli- Villars regularization and discrete light-cone quantization 
to the nonperturbative solution of a (3 + l)-dimensional model held theory. 
The matrix eigenvalue problem is solved for the lowest-mass state with use 
of the complex symmetric Lanczos algorithm. This permits the calculation of 
each Fock-sector wave function, and from these we obtain values for various 
quantities, such as average multiplicities and average momenta of constituents, 
structure functions, and a form factor slope. 
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I. INTRODUCTION 



One of the most challenging problems in particle physics is the computation of the 
spectrum and physical properties of bound states in quantum field theory. The main tool 
presently used for such nonperturbative computations in quantum chromodynamics is lattice 
gauge theory [Q, which has been highly successful for determining hadron spectra. However, 
the computation of dynamical properties, such as CP violation in weak transition matrix 
elements or the shape of the distributions measured in deep inelastic scattering is difficult 
using standard lattice methods. 

Light-cone Hamiltonian diagonalization methods |Q appear to provide a number of at- 
tractive advantages for solving nonperturbative problems in quantum field theory, including 
a Minkowski space description, boost invariance, no fermion-doubling, and a consistent Fock 
state expansion well matched to physical problems in QCD; however, thus far, full dynam- 
ical solutions based on light-cone Hamiltonian diagonalization have been primarily limited 
to one-space/one-time models. One promising approach is the transverse lattice which com- 
bines light-cone methods in the longitudinal light-cone direction with a spacetime lattice for 
the transverse dimensions. 

In recent work ^ we have shown that a model field theory in 3 + 1 dimensions can 
be solved using discrete light-cone quantization (DLCQ) a light-cone Hamiltonian 

diagonalization method, together with Pauli-Villars regulation of the ultraviolet The 
particular model theory which we constructed has an exact analytic solution by which the 
DLCQ results could be checked, for both accuracy and rapidity of convergence. The model 
was regulated in the ultraviolet by a single Pauli-Villars boson, which was included in the 
DLCQ Fock basis in the same way as the "physical" particles of the theory. The two bare 
parameters of the model were then determined by fits of observables to chosen values. 

Here we shall extend this combination of DLCQ and Pauli-Villars regularization to a 
more realistic model which mimics many features of a full quantum field theory. Unlike the 
analytic model which contained a static source, the light-cone energies of the particles in the 
new model have the correct longitudinal and transverse momentum dependence. Although 
an analytic solution of the new model is no longer available, the numerical convergence 
of the discretized light-cone solutions is found to be quite rapid, and the structure of the 
solution for the lowest-mass eigenstate is readily obtained. In particular, we can calculate the 
light-cone wavefunction of each Fock-sector component, and from these we can compute the 
values for various physical quantities, such as average multiplicities and average momenta 
of constituents, bosonic and fermionic structure functions, and a form factor slope. 

A distinct advantage of our approach is that almost all counterterms are generated 
automatically by the Pauli-Villars particles and their imaginary couplings. This can be 
explicitly checked for consistency in perturbation theory. For nonperturbative calculations 
we conjecture that the same number of Pauli-Villars fields will be sufficient to regulate the 
theory. This does appear to be the case here and in the work reported previously 0. An 
alternative procedure has been proposed and explored by Wilson, Perry and collaborators 
P]; they use a similarity transformation to generate effective Hamiltonians perturbatively 
which can then be diagonalized in the valence Fock sector. 

In our approach one can obtain the full set of Fock-sector wave functions for the lowest- 
mass eigenstate. This contrasts with other DLCQ calculations in 3 + 1 dimensions [^-[TT 
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where the number of particles was severely limited from the outset and effects of higher Fock 
sectors can only be estimated. The DLCQ calculation by Wivoda and Hiller though 
untruncated, did not construct counterterms in a way that can be systematically extended 



to other theories. In our case, a Tamm-Dancoff truncation in particle number can be 
applied, and the impact of the truncation can be studied and understood. 
Our notation is such that we define light-cone coordinates |]I4| by 

x^ = x^ + x^, x^ = {x^,x^). (1.1) 

The time coordinate is taken to be x^ . The dot product of two four- vectors is 

p ■ X = -{p^x^ + p^x ) — p± ■ x_i . (1.2) 

Thus the momentum component conjugate to x~ is p'^, and the light-cone energy is p~ . We 
use underscores to identify light-cone three-vectors, such as 

p = (p+,px). (1.3) 

For additional details, see Appendix A of Ref. |0 or a review paper |Q. 

The model which we study is defined in Sec. ||. There we also list and define various 
quantities which we will compute from the eigensolution, including structure functions and 
distribution amplitudes, average multiplicities, and average momenta. The numerical meth- 



ods, including the DLCQ procedure, and the results are discussed in Sec. |T|. Section 
contains some concluding remarks and plans for future work. 



II. A MODEL WITH A DYNAMICAL SOURCE 

We shall consider a field-theoretic model where one particle, which we take to be a 
fermion of mass M, acts as a dynamical source and sink for bosons of mass /i. The model is 
only slightly more complicated than the analytically soluble model considered in Ref. |^ , the 
key difference being that here the fermion has a proper, momentum-dependent light-cone 
energy. Another difference is that the vertices do not include the momentum ratios which 
were introduced in |Q to control end-point behavior; the restoration of fermion dynamics 
makes such factors unnecessary. The theory is still regulated by a single Pauli-Villars boson 
with imaginary couplings^ and a mass fii. The light-cone Hamiltonian (or mass-squared 
operator) i^LC = P^P^ — P± is, in the P± = frame, 

_ f dp+<fp^ ( m'+pI p+ \ V , t , 



One could use an Hermitian form and negative metric to implement Pauli-Villars regularization, 
but the complex symmetric form is what is known to work well with the numerical method we 
have chosen. 
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+ 



dq'^(Pq± 
167r'^g+ 



/• dptd^p±i r dptd^p±2 f dq+d'^q^ ^ x 



1Qti^P2 

where aj, and aj^ are creation operators for the fermion source, the physical boson, and 
the Pauh-Villars boson, respectively. The operators obey the usual commutation relations 

{V,&p'^'} = 167r^]9+5(p-p')'^<T<7' , 
167r3g+(5(g - q) , 
167r3g+5(g - q) 



t 



ai„,a 



(2.2) 



The MqP^/P^ counterterm is inserted to cancel a logarithmic dependence on the Pauli- 
Villars mass which arises from the one-loop self-energy integral 




d'q± 



P+/P+ 



M2+(px+qx)2 

{p+-q+)/P+ 



Q+/P+ 



— P-V term 



ln(/ii//i). (2.3) 



This model Hamiltonian is distantly related to the Yukawa Hamiltonian |T5|, to which one 
might also eventually apply the techniques used here. 

The bare parameters g and Mq are to be fixed by fitting physical properties of the lowest 
massive eigenstate. This is a dressed fermion state which we write as 



n,n\ 



dp'^d^p± f dq^d^q^. 



n 



i=l 



n 



dr^d^r^j 



(2.4) 



X 



5{P~p-Y,gi-Jl': 



/ (n,ni) 



and normalize according to 

<|.;t . $^ = 167r^P+5(P' - P) . 
The individual amplitudes must then satisfy 



(2.5) 



E n y dqtd'q^, n J dr^d^^j 



n,ni I 



(--)(g^,r,;P-E^.-Ei:. 



(2.6) 



The eigenvalue problem is 



(2.7) 
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This is equivalent to the following coupled set of integral equations for the amplitudes: 



X 



E 



(n,ni) 



(2.8) 



E 



1 



(n— l,ni ) 



+Wni + 1 y 



{q,,rj,r,p-r) 



+ 



E 



(n.ni — 1) 



(^j) Hi; • • • jEj-l; Ej + l; • • • , Lni^ P + Hj) 



with X = p'^ /P~^, Hi = /P^, and Zj = r^/P"*". 

For fixed M, the eigenvalue problem itself is a condition on the bare parameters. A 
convenient choice for the second condition is the value of an expectation value involving the 
boson field (j){x); we use (:0^(O):) = $J. :0^(O): $o-, which corresponds to the expectation value 
for the sum of 2/?/j for physical bosons. For the soluble model in Ref. it was shown to be 
closely tied to the coupling g, as can be seen in Eq. (3.11) of that paper. Most importantly, 
it can be computed rather quickly from a sum similar to the normalization sum 



„ "1 „ / n 

(:0'(O):) = En=i,n,=on? J dqtd'q^^H J dr+d\^, ^ 



.=1 <lt/P' 



(2.9) 



These two conditions are sufficient to fix g and Mq. 

With the two parameters of the model now fully determined, we can compute other 
quantities as predictions. These are all obtained from the primary output, which is the 
set of wave functions for the different Fock sectors. We will compute the slope of 

the no-flip form factor of the fermion, structure functions for bosons and the fermion, the 
distribution amplitude for the physical boson, average momenta, average multiplicities, and 
a quantity sensitive to boson correlations. The form factor slope -F'(O) is given by [^] 



F'W = E ny dqtd\^,Jl j dr^dh^ 



(2.10) 



X 



T.. f vL + E fvi, 



(n,ni) 



X( 



(n,ni) 



(g,,£j;£-E^i-Eij)- 



A related form. 
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(2.11) 



n.n\ I 



E fVx.0("-ng„i:,;£-E^.-Ei:.; 



+ E |v^.'/>^--Hg^,r,;P-Ei,-E^ 



is better computationally. It is obtained from (|2.10|) via integration by parts. If a momentum 
cutoff is present, there are surface terms, but these will vanish at infinite cutoff. 
The physical boson structure function is defined as 



fsiy) = En,n, m I dqfd'qu Yif I dr^d^^, ^ 6{y - qt / P' 



(2.12) 



1=1 



X 



{n,ni) 



(g,,i,;Z-E^i-E^i 



The fermion and Pauli-Villars structure functions fpix) and fpy{z) are defined analogously. 
The normalization of each is such that the integral yields the average multiplicity 



(2.13) 



{ns) = fBiy)dy, {ripv) = fpv{,z)dz . 
Jo Jo 

The average momentum carried by each type is also given by an integral 

{y)= f\fB{y)dy, {z)= (\fpv{z)dz. (2.14) 
Jo Jo 

As a measure of the correlations in the multiple-boson Fock sectors, we compute the covari- 
ance (i/i?/2)n>2 - {y)l>2 where 



(2.15) 



X 



/.(n.rii) 



\%,r-,P-Y^q.-Y.rj 



and {y)n>2 is the same as (y) except that only states with two or more bosons are included. 
We also compute the distribution amplitude [|1^ given by (f{y) = J d'^q±(p^^'°\y, q±). 



III. NUMERICAL METHODS AND RESULTS 

A. Discretization and diagonalization 

We discretize the coupled integral equations and the formulas for quantities such as the 
form factor slope in the standard DLCQ manner Integrals are approximated by discrete 
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sums and derivatives by finite-differences. Because of tlie Pauli-Villars regulation, tlie tlieory 
is ultraviolet finite. However, in order to have a finite matrix problem, we limit the range 
of transverse momentum by imposing a cutoff on each constituent's invariant mass 

!li±^<A', (3.1) 

where is the physical mass of the constituent. (Later, we study the large A limit.) The 
longitudinal momentum, always being positive, has a natural finite range. 

Given the length scales L and L±, the discrete momentum values are taken to be 

p^-*-n, p± ^ {—n.^,—ny) , (3.2) 

Li 1j± h±_ 

with n even for bosons and odd for fermions. The differing values of n correspond to use of 
periodic and antiperiodic boundary conditions, respectively, in a light-cone coordinate box 

— L < x^ < L , —L± < x,y < L± . (3.3) 

The total longitudinal momentum is used to define an integer resolution [^] = ^P^- 
The positivity of the longitudinal integers n implies that the number of particles in any 
Fock sector is limited to ~i^/2. The integers and Uy range between limits associated 
with some maximum integer N± fixed by L± and the cutoff A, such that N±it/L± is the 
largest transverse momentum allowed by the cutoff. 

The integral equations and other physical objects are independent of L, a feature of 
boost-invariance in DLCQ. The limit L — > oo is replaced by the limit K —>■ oo. The 
momentum-space continuum limit is reached when both K and N± become infinite. The 
momentum-space volume limit A^ ^ oo is taken after the continuum limit. 

Weighting factors are included in the sums that approximate integrals in order to incor- 
porate boundary effects induced by the invariant-mass cutoff. For a discussion of how these 
factors are constructed and used, see Ref. p. 

Typical basis sizes are given in Table |I[ The present calculations, which use a single four- 
processor node of an IBM SP, are limited to ~11 million states. The Hamiltonian matrix is 
extremely sparse, so that the lowest-mass state can be efficiently extracted with use of the 
Lanczos algorithm for complex symmetric matrices ||18|J5[|. The analytic solution for the 



soluble model discussed in Ref. is used as an initial guess for the Lanczos procedure. 

Before invoking the Lanczos algorithm, the eigenvalue problem is rearranged so that 
— 1/g is the eigenvalue. This allows computation of g given a fixed value for M and a guess 



for Mq. The iterative Brent-Miiller algorithm is then used to find the value of Mq that 



brings (:0^(O):) into agreement with its chosen value. 



B. Results 



Most of the calculations reported here use the parameter values = /x^ and (:0^(O):) = 
1. These choices correspond to a relativistic, weak-coupling regime. Because of the weak 
coupling, the number of Fock sectors can be truncated to include no more than four bosons 
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TABLE I. Basis sizes for DLCQ calculations with parameters = fi'^, fi^ = and 
= 50 fj?. The numbers of physical states are in parentheses. 











K 






9 


11 


13 15 17 


5 




54100 


95176 


386140 1553 576 6 816 394 






(28 065) 


(66 371) 


(232 400) (1038 070) (4 972 065) 


6 




126 748 


536 758 


2 907158 4 935 510 






(69 245) 


(391511) 


(2107 688) (3 013 689) 


7 




519 325 


1317 392 


10 080 748 






(276 299) 


(1 008 539) 


(7 272 134) 
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1 165 832 


5162 002 








(687 394) 


(4 140 491) 
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2 268 535 










(1437 647) 
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5 850 335 










(3 585 752) 






TABLE II. Fock sector probabilities 


J ^(n,?ii) 2 j-jra -j-jrai ^hgrg n ig the number of 


physical bosons and rii the number of Pauli-Villars bosons. The numerical and physical parameters 


are 


K = 17, N± = 7, M2 = 


= /i^, fij = Wfi 


2, A2 = 50;u2, and (:0^(O):) = 1. The total number of 


bosons n + ni is limited to 


a maximum of 4. Probabilities smaller than ~ 10 ^ are not resolved 


with 


any accuracy. 






n\ni 





1 


2 3 4 







0.8515 


0.0115 


0.8-10-^ ~ 10-^° ~ 10"^*^ 


1 




0.1333 


0.0005 


~ 10-'^ ~ 10-^2 


2 




0.0036 


0.4-10"^ 


~ 10-10 


3 




0.3-10-^ 


-10-8 




4 




~ lO"'^ 







without any discernible effect, as can be seen from the Fock-sector probabihties hsted in 
Table most weak-coupling calculations were done with this truncation in order to increase 
the available momentum resolution. For comparison, we have also done some study of other 
regimes. 

Table shows values of various quantities, extrapolated from longitudinal resolutions 
-fC = 9 to 19 (or even 21) and transverse resolutions A^_l = 5 to 10 for small K and to 6 or 
7 for large K. These include the bare coupling g, the renormalization parameter Mq, the 
bare fermion probability l^/'oP? the slope of the form factor F'{0), the average multiplicity 
(ub), and a parameterization of the structure function fsiy) = Ay°-{1 — y)^ (which is an 
excellent fit). Each is shown as a function of the cutoff and the Pauli-Villars mass /ii. 
The extrapolations were done by fitting to the form a + (3 / K"^ + 'y/N]_] most quantities are 
slowly varying with respect to resolution. The range of values obtained for -F'(O) correspond 
to a dressed-fermion radius J—QF'{0) on the order of 0.2 



8 



TABLE III. Extrapolated bare parameters and observables. The physical parameter values 
were = /j.'^ for the fermion mass and (:(/)^(0):) = 1. 





/if = 5/i2 


/^? = 


10/x2 




20/i2 


(A//.)^ 


12.5 


25 


50 


25 


50 


50 


100 




21.4 


17.7 


16.3 


17.8 


16.0 


16.0 


15.5 




1.26 


1.10 


1.10 


1.48 


1.4 


1.8 


1.9 




0.82 


0.83 


0.84 


0.85 


0.86 


0.87 


0.87 


1 nn;/2 ^Ur\\ 

— lUu/i X* \\J j 


1.04 


0.78 


0.66 


n 79 


u.oy 


u.oy 


U.Oi 


{iT'B) 


0.18 


0.15 


0.14 


0.15 


0.14 


0.13 


0.13 


(y) 


0.077 


0.062 


0.057 


0.062 


0.056 


0.056 


0.053 


{yiy2)n>2 - {y)l>2 


1.1 • 10-3 


6 • 10-4 


6 • 10-4 


6 • 10-4 


6 • 10-4 


6 • 10-4 


5 • 10-4 


A 


9.39 


4.21 


3.00 


4.15 


2.77 


2.7 


2.4 


a 


1.90 


1.50 


1.36 


1.48 


1.31 


1.29 


1.26 


b 


2.95 


2.54 


2.32 


2.53 


2.26 


2.24 


2.14 



The table shows that the renormalization parameter Mq is the only quantity strongly 
dependent on the Pauli-Villars mass. This is to be expected because of its role in the 
self-energy counterterm. One might argue that F'{0) is also strongly dependent; however, 
any apparent variation with /if is largely due to differences in cutoff values and transverse 
resolution. Although -F'(O) will ultimately become independent of and A^_l, it is sensitive 
to these in the ranges where we calculate. The table also shows that the estimate of J^iVi << 
1 by Glazek and Perry is justified, in that the expectation value (y) is found to be small. 

A sample boson structure function is plotted in Fig. |1]. The figure also shows how 
well the form 74?/"(l — y)^ fits the numerical results and how insensitive Jb is to numerical 
resolution, something which was also observed for the model considered in Ref. ||^. The 
transverse and longitudinal dependence of a two-body amphtude are shown in Figs. ||, |^ 
and 1^. A particular transverse cross section of the two-body amplitude is presented in 
Fig. ^; these results correspond to fixed values of the transverse scale L±_ and are remarkably 
consistent. Figure ^ shows the dependence of the boson structure function. A fermion 
structure function and a Pauli-Villars boson structure function are plotted in Figs. |^ and 
^. The parameter values are the same for both. The skewing of the Pauli-Villars particle 
momentum distributions to high longitudinal momentum fractions reflects the heavy mass 
of the Pauli-Villars bosons. 

Other values for the physical parameters M and (:(/)^(0):) have also been considered. A 
summary of extrapolated quantities is given in Table The associated structure functions 
fsiy) are shown in Figs. |^ through Distribution amplitudes are displayed in Figs. ^ 
and 0. For values of M larger than /i we have found the form A?/'*(l — y)^e~'^y to allow 
a noticeably better fit to fsiy)- For (:0^(O):) = 5 the maximum number of bosons was 
increased to 5. The numerical resolutions ranged from 9 to 21 for K and from 5 to as much 
as 10 for N±. 

The extent to which the fermion source is dressed by the bosons is directly determined 
by the mass ratio M/ fi and the coupling strength. The latter is tightly correlated with the 
chosen observable (:0^(O):). As the ratio M//i is tuned, the boson structure function fsiy) 
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0.30 




FIG. 1. The boson structure function /b at various numerical resolutions, with M = ji, 
{■.(f>'^{0):) = 1, = and /xf = The solid line is the parameterized fit, Ay"'{l — y)^. 
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TABLE IV. Same as Table 



III. 



, but for different or (:(/)^(0):) values. 





{■4\0):) = 1 




{M/^,f 


0.1 


5 


10 


1 




10 


10 


10 


10 




50 


100 


100 


50 




15.1 


18.1 


19.0 


44.5 




1.39 


1.66 


1.60 


10.1 




0.83 


0.89 


0.90 


0.41 


-100/i2F'(0) 


2.0 


0.14 


0.07 


6.7 


{^b) 


0.16 


0.10 


0.09 


0.62 


(y) 


0.073 


0.032 


0.024 


0.24 


(yiy2>n>2 - {y)l>2 


7 • 10-4 


3 • 10-4 


3 • 10-4 


8 • 10-3 


A 


1.0333 


5.2548 


7.5519 


9.0847 


a 


1.0512 


1.3191 


1.3339 


1.0256 


b 


0.8678 


2.5430 


1.7151 


2.1580 


c 





2.2730 


4.9870 






shifts dramatically. A relatively small boson mass shifts the peak in fsiv) to small boson 
momentum fractions, as shown in Fig. |T^. A large mass shifts the peak to central values 
of y and significantly raises the constituent density at large y, as illustrated in Fig. |^. An 
increase in (:(/)^(0):) increases the coupling and increases the probability for a large number 
of constituents. Analogous changes occur for the distribution amplitude. Comparison of 
Tables |ITl| and fV\ shows that the average number (ns) increases significantly when (:0^(O):) 
is changed from 1 to 5. 



IV. CONCLUSION 

We have successfully computed the Fock-sector wave functions which fully describe the 
lowest-mass eigenstate of a field-theoretic model Hamiltonian (|2.1| ) in physical three space 
and one time dimensions. From these wave functions we have extracted several interest- 
ing quantities to show that numerical convergence is under control and that Pauli-Villars 
regularization leads to sensible results. The size of the momentum-state basis required is 
large but manageable for present-day computing machines. Larger bases could be used by 
expanding to more than one node, although one then pays the price of message-passing 
overhead. 

For the model discussed here there are still interesting calculations which might be done. 
One could look at excited states in the one-fermion sector that we have explored, or consider 
other sectors, such as the two-fermion sector. Extension to two flavors, particularly with 
very different masses, should yield some understanding of light systems with heavy intrinsic 



constituents, which could have some relevance for intrinsic charm ^1 . 

Beyond this model there are, of course, many possibilities. A solution of Yukawa the- 
ory P2[, in a no-pair approximation or eventually in full, would be the most immediate 



nontrivial extension. Applications to quantum electrodynamics, to positronium |10] or the 
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electron's anomalous moment |Tl| in particular, would be quite natural. Direct applica- 



tion to quantum cliromodynamics (QCD) may be problematic; however, a supersymmetric 
conformally-invariant form of QCD could lend itself to the spirit of the approach, in that 
heavy superpartners in a broken supersymmetry should provide the needed ultraviolet can- 
cellations. 
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FIG. 2. The one-boson amplitude (p^^''^^ as a function of longitudinal momentum fraction y and 
one transverse momentum component Qx in the Qy = plane. The parameter values are K = 21, 
= 7, fij = 10^^, = 25/x2, and (:(/'2(0):) = 1. 
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FIG. 3. The boson-fermion two-body amplitude at zero transverse momentum, with K = 21, 
N± = 7, (:^^(0):) = 1, = and /xf = 10/i^. The normahzation is arbitrary. 
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FIG. 4. Cross sections of the boson-fcrmion two-body amplitude taken at varying longitudinal 
momenta and at fixed qy = 0, with K = 21, iVj_ = 7, (:0^(O):) = 1, = and /xf = 10//^. 

The peaks are normalized to be equal at Qx = 0. 
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FIG. 5. A cross section of the boson-fermion two-body amplitude taken at fixed longitudinal 
momentum fraction y = 4/9 and at fixed Qy = 0, witfi ii' = 9, {: (j? {Q) :) = 1, and /if = 
The cutoff and the transverse resolution Nj^ are varied to keep the transverse scale L±_ fixed 
at one of the following values: 1^ (black), \/2^ (gray), and 2^ (white). Different symbol shapes 
correspond to different values of A^j^. The peaks are normalized to be equal at Qx = 0. The points 
at zero amplitude mark the transverse range, which is set by the cutoff. 
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FIG. 6. The boson structure function /s(2/,q±) with K = 21, = 7, (:0^(O):) = 1, 
= 25;u^, and /xf = The transverse momentum is varied with Qy fixed at zero. 
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FIG. 7. The fermion structure function fp with K = 21, N± = 5 to 7, (:0^(O) :) = 1, 
= and /xf = Each N± value yields essentially the same result. 



20 




0.000 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 8. Same as Fig. ^ but for the Pauli-Villars boson structure function fpv- 
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FIG. 9. Same as Fig. |l|, but for = 0.1^^ 
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FIG. 13. Comparison of distribution amplitudes (p{y) = / (fq±4>^^''^\y,q±)- Various values 
are considered for the fermion mass M, with {:(jp{0):) = 1. The values of the numerical parameters 
are K = 19, N± = 5, = 50/^^ (except for = 5fi^ and lO^t/^ when = lOO^u^), and fij = Wfi"^. 
The lines simply connect the computed points, to guide the eye, and the absolute normalization is 
arbitrary. 
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